Orbital liquid in three dimensional Mott insulator: LaTiO^ 
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We present a theory of spin and orbital states in Mott insulator LaTiOs . The spin-orbital superex- 
change interaction between d^{t2g) ions in cubic crystal suffers from a pathological degeneracy of 
orbital states at classical level. Quantum effects remove this degeneracy and result in the formation 
of the coherent ground state, in which the orbital moment of t2g level is fully quenched. We find a 
finite gap for orbital excitations. Such a disordered state of local degrees of freedom on unfrustrated, 
simple cubic lattice is highly unusual. Orbital liquid state naturally explains observed anomalies of 
LaTiOs- 
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Experiments on transition metal oxides continue to 
challenge the theory - Keimer et al. have recently re- 
ported on dynamical quenching of t2g orbital angular mo- 
ments in Mott insulator LaTiO^ which is not associated 
by any detectable orbital/ Jahn- Teller ordering. More 
specifically, magnon spectra are found to fit a spin one- 
half, nearly isotropic Heisenberg model on cubic lattice. 
At the same time, anomalous x-ray scattering study has 
not revealed any evidence of static orbital order. These 
observations indicate that t2g orbital degrees of freedom 
of Ti'^'^ ions are mysteriously missing in a low energy 
states of LaTiO^. Keimer et al. suggested a picture of 
strongly fluctuating orbitals to reconcile their results. 

Dynamical quenching of local orbital moments in a 
periodic, simple cubic lattice formed by Ti^+ ions in 
LaTiOs is fascinating. This poses a serious problem for 
the canonical Goodenough-Kanamori picture of succes- 
sive orbital and magnetic orderings - the guiding idea 
in the modern theory of orbitally degenerate transition 
metal oxides (see for review [^,^). Loosely speaking, a 
lifting of the orbital degeneracy without symmetry break- 
ing is formally similar to the spin moment quenching in 
periodic Kondo lattices. The novelty of LaTiO^ is how- 
ever that it is an insulator. At the same time, the ideas 
of quantum disorder due to the weak connectivity or ge- 
ometrical frustration (like in pyrochlore systems, see 
for instance) do not apply here, either. 

As discussed below, there are several physical reasons 
for the orbital disorder in LaTiO^, and the key factor is 
a special symmetry of t2g -superexchange interaction on 
cubic lattice. This results in a quantum resonance be- 
tween orbital levels, which removes degeneracy present 
at classical level by the formation of resonating bonds in 
orbital sector - a realization of Anderson's resonating va- 
lence bonds (RVB) idea H] in a three dimensional (3D) 
insulator with help of orbital degrees of freedom. As for 
spin moments, they show a long-range order as expected 
on general grounds. However, a staggered moment is con- 
siderably reduced, since short- wavelength magnons are 
actively involved in the local resonance of exchange bonds 
- in fact, it is the composite spin-orbital excitation which 



plays a crucial role in the theory presented. 

We begin with discussion of the t2g -superexchange 
interaction. In terms of fermions ai,cr, bi^a, Ci^a corre- 
sponding to t2g levels of yz, zx, xy symmetry, respec- 
tively, the hopping term in the Hubbard model reads 

as —t {a\^„aj,a + bl^^bj^a- + H.c.^ for the bonds along c- 

direction in a cubic crystal. Correlation energies in dou- 
bly occupied virtual states of Ti^^ ion are specified as 
follows: U for electrons on the same orbital; U + Jh, 
if they occupy different orbitals and form a spin-singlet; 
and U — Jh for a spin-triplet charge excitation. The 
relation U = U — 2Jh holds in the atomic limit. As 
Jh << U usually, the dominant part of the superex- 
change interaction in a cubic lattice is: 
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On c-axis bonds the orbital structure of the exchange 
"constant" is given by (hereafter, the energy scale At'^/U 
is used): 



a%iblai + b]ajalb 



(2) 

(b) 



Similar expressions can be obtained for J^"''^ and 
by replacing orbital fermions a, b in this equation by b, c 
and c, a pairs, respectively. We notice, that in Eq.(2) 
and below fermions have only orbital quantum number 
(a,5, c orbitons), since a spin component of the original 
fermions is represented in Eq.(l) by spin one-half oper- 
ators, as usual. We shall focus now on the Hamiltonian 
(1), and discuss later the effects of finite Hund coupling 
corrections (of the order of Jh/U). 



As evident from Eq.(l) (consider (S^Sj 



-1/4), the 



classical Neel state is infinitely degenerate in orbital sec- 
tor, refiecting that cubic symmetry is respected by the 
Neel state. This emphasizes a crucial importance of 
quantum effects, as has been pointed out first in the con- 
The bond directional geometry 



text of gg-systems ] 
of Cg orbitals offers a solution 



in that case: eg or- 



bital frustration is resolved by order from disorder mech- 
anism selecting a particular (directional) orbital configu- 
ration, which maximizes the energy gain from quantum 
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FIG. 1. (i) Spin-orbital resonance between two (left and 
right) states of the c-bond pair, (m) Self-consistent Born ap- 
proximation for the fermionic orbital self-energy, (in) Di- 
agram contributing to the bond order parameter x- Solid 
(wavy) lines represent orbitons (magnons). 

spin fluctuations. t2g orbitals are however not bond ori- 
ented (they are all planar), and not much spin fluctuation 
energy can be gained by any pattern of static orbital or- 
derings. The solution of the <2g-problem proposed here 
is different: that is an idea of S'C/(4) spin-orbital excita- 
tions recently discussed in a context of one dimensional 
models 

Reflecting the geometry of t2g orbitals, every bond in 
Eq.(l) is represented by two equivalent orbitals (such a 
symmetry is absent for Ising-like eg orbitals). This is 
also evident from the representation of Eq.(2) in terms 
of orbital pseudospin t^^ and density nj^^ — riia + nib 
operators acting on the (a, 6)-doublet subspace: J-^j^^ = 

2 {^IbTib + \<bKi) ■ Suppose that nj,^ = 1; one then 
clearly observes an orbital SU (2) symmetry of exchange 
"integrals". As we know from spin 5t/(2)x orbital SU{2) 
models ^|-|lo| , the exchange energy is gained in that case 
due to the resonance between degenerate local configu- 
rations (spin-singlet x orbital triplet, and spin-triplet x 
orbital singlet), and elementary excitations are a mixed 
spin-orbital S'[/(4) modes. Fig.li'ij illustrates the idea of 
correlated spin-orbital fluctuation. In the present prob- 
lem, spins and orbitals are different: spin sector is ef- 
fectively half-flUed (f^i,| + = 1) and must show a 
long-range order in 3D; on the other hand, the orbital 
occupation number n^^ is not conserved but rather fluc- 
tuates around the average value 2/3. Because of the 
presence of third orbital, one has a less severe constraint 
n-i^a+ni^b + ni^c = 1 in orbital sector. Even though a true 
SU{A) excitation cannot develop here because of spin- 
orbital asymmetry, the analogy with ID models suggests 
the way how to optimize the exchange energy: That is by 
the formation of a virtual 5t/(4) resonance which can be 
viewed in the present context as a local excitation com- 
posed of magnon and orbital fluctuation. This resonance 
removes the orbital frustration dynamically, and the dis- 



ordered state of orbitals is precisely what is required to 
amplify this mechanism in all three directions of cubic 
lattice. 

Technically, the main ingredient of the theory is the 
three-particle bond variable, Xa'^^ — a|aj(s| + Sj)- Here, 
a denotes an appropriate fermionic orbiton, say a or 6 or- 
bitons for c -bond pair {ij), while being a magnon exci- 
tation about the Neel state. This is an analogy to SU (4) 
flavor bond variable |^ in a spin ordered state. The bare 
magnon dispersion is simply given by the Heisenberg in- 
teraction with an average exchange constant J = {J^j'^), 
which will be calculated below. On the contrary, the 
orbital dynamics is exclusively due to the three-particle 
resonance, since spin-only prefactor (S^Sj -I- 1/4) is al- 
most zero —0.05) in 3D. We may present Eq.(l) as 
^^SE = ^sp + ^int + const., whcrc Hsp describes bare 
magnons, and Hint = -\ Xa'^^X^^*'' • The strat- 
egy is then to decouple _ffi„f in terms of bond vari- 
ables. The solution with uniform bond amplitude % gives 
H-mx. ^ Ha + Hb + He + const., where 

Ha = -X-Y. (4 + -'i) + H-C-} ■ (3) 

{ij)ba 

Here the summation is restricted to (6, c)-planes. Similar 
expressions hold for b and c orbitons. Since the filling 
factor is only 1/3 per orbital flavor, the orbital sector is 
far from nesting conditions, thus a mean-field uniform 
solution in the spirit of large N -theories [0 is believed 
to be a reasonable starting point to describe a disordered 
state of orbitals. Eq.(3) describes an effective hopping 
of orbitons accompanied by simultaneous magnon exci- 
tations. Thus, we have mapped the superexchange in an 
undoped system to an effective x — J model, which in a 
momentum space reads as 

^x-j = J2J2 ^-4? («k+p«k5p + H.c.) + J2 ^p4^p- 

kp a p 

(4) 

Here M^^ = — 4x(7^"^pUp -I- ll^^Vp), while form- factors 

are 7^^^^ = {cy + c.)/2, 7^ = + c,)/2, 7^'^ = 
(c^; + Cy)/2, where cx — cosk\. Bare magnon disper- 
sion LOp = 3J(I — 7p)^^^ with 7p = [cx + Cy + Cz)/^^ and 
Mp,fp are conventional Bogoliubov transformation coef- 
ficients for magnons in the cubic lattice. It is noticed 
that the coupling constant in Eq.(4) vanishes at p = 
limit (spin conservation is respected), and it is the short- 
wavelength magnons which are important. The model is 
similar to the t — J model for the doped Neel state of 
cuprates (see [|l2|,^ for comparison) , yet the parameters 
X and J have to be self-consistently determined. It is 
noticed, that there are three branches of 2D fermions, 
and the dimensionality of magnon sector is three in the 
present case. 
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We calculate first the fermionic spectrum within a self- 
consistent Born approximation. Different from holes in 
the t — J model, an orbiton creates a magnon on the 
site it arrives at, and eliminates a magnon while leav- 
ing the site (see Eq.(3)). Therefore orbiton motion con- 
tains a coherent component even in the Ising limit for 
magnons. The latter approximation is used to simplify a 
momentum integration in the self-energy, Fig.l(M). One 
then obtains i?eEif''(k) = — '«u;(7k"')^- Here k,^ — 
(4x)2/o°°dep(0/(^ + 3J-C.), where p{i) = (pk(0>k 
is fermionic density of states (DOS), has similar 
structure: = [Axf d^p{^)/{uj -f 3 J - where 

p(0 — (Pk(0(7k"'')^)k- One observes a two-sublattice 
structure of the orbiton dispersion (imposed by spin or- 
der) as expected from analogy with the t — J model 
p^ . We estimate the fermionic (unrenormalized) mass 
as m ~ I/kq, and density of states as p(0) ~ I/ttko. 
Each orbiton forms its own 2D fermi-surface (FS). Say, 
for a orbitons the FS consists of two almost circles around 
(0, 0) and (tt, tt) points in a {ky, kz) plane. Assuming con- 
stant DOS within the interval of width Worb = l/p(0), 
we find kq — ^xl\f^i ^^nd Worb — '^y^X- Worb is the 
energy scale of orbital fluctuations jlj]. Further, we 

estimate the bond amplitude x = {'A'^ji^l + ^i))- 
the Ising limit for intermediate magnons, Fig.l(m) gives 
X ~ (4/3V^)((7^''^)2)j.s- For = 1/3 we estimate 

{{i^''Y)fs — 0.3. Spin stiffness in the present model is 
controlled by J. The mean-field value 2/9 for J follows 
from Eq.(2).This number is reduced to J ~ 0.16 due to 
the corrections shown in Fig.2(i). The physics behind 
this reduction is that SU{A) resonance induces some fer- 
romagnetic component in spin interactions. 

Thus, we have fixed basic energy scales (in units of 
At^lU): Worb - 1-6 for orbital sector, and 3J ~ 0.48 
for magnon bandwidth. The spin stiffness is relatively 
small because Pauli principle is relevant only if electrons 
occupy the same orbital, and that probability is reduced 
at large orbital degeneracy. Main energy gain stems from 
magnon-orbital resonance: Emt — —Qx^ — —0.31. (Spin- 
only contribution is negligible: Esw — —0.02). This 
number should be compared with energy gain that one 
has with statically ordered planar orbitals: ~ —0.17. 

We discuss now physical consequences of the theory in 
context of experimental data on LaTiO^ [|l|. Observed 
magnon dispersions of cubic symmetry are very natu- 
ral in above picture of orbitals fluctuating faster than 
magnons. The energy scale At^ /U ~ 100 meV, required 
to flt Jexp — 15.5 ± ImeV, is also reasonable in view 
of values f7 - 4eV |l|, and i - 0.3eV |l|]. This 
scale means that orbital fluctuation energy in LaTiO^ 
is Worb ~ leOmeV. 

Keimer et al. observed very small anisotropy gap, indi- 
cating that effects of conventional relativistic spin-orbit 
coupling are strongly suppressed. This is in fact an un- 
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FIG. 2. (i) Interaction corrections to the spin exchange 
constant J. (ii) Magnon self-energies. 



avoidable consequence of our theory, since the ground 
state has no orbital degeneracy. We may calculate spec- 
tral density of orbital angular momentum fluctuations. 
It is given by the fermionic orbiton excitations ( "orbiton 
Stoner continuum"). Spectral density of the local sus- 
ceptibility of angular moments If — i{a% — a)i van- 
ishes in the static, a; = limit: S{uj) — 2/9^(0)^. This 
implies complete quenching of orbital moments at low 
energies. Making parallel with Jahn- Teller (JT) impu- 
rity physics jl^, we may say that "effective Ham reduc- 

r - 1/2 

tion factor" , C,{Cj) ^ \ S{uj)duj > ^ ujjWorb is fre- 
quency dependent, and the angular moments disappear 
from low-energy physics linearly in energy. Also, the ratio 
^sojWorb ^ 10"^ (Ago ^ 20 meV [0) means that rela- 
tivistic spin-orbit coupling induced corrections are small; 
we estimate the g-value shift Ag/g 0.06, consistent 
with 0. 

Despite three dimensionality of magnon spectra, the 
staggered moment in LaTiOs is small, i.e. 0.45/1 _b - 
an obvious problem for a spin-wave picture. The present 
theory resolves this difficulty: The intensity of spin Bragg 
peak is partially taken away by a quantum magnon- 
orbiton resonance and redistributed over finite frequency 
region. In other words, fluctuating orbitals generate ad- 
ditional quantum spin fluctuations in the ground state. 
We have calculated a spin moment reduction in a sim- 
ilar way as it was previously done in the t — J model 
p8[ . Accounting for lowest order interaction corrections 
(Fig.2(M)) to the magnons, one finds: 



2 f-' 1^ (wp +4k' - Ck)-^ 



kk' 



Skk' 



i^pi^p + Ck' - Ck) 



(5) 



where ^k is the a orbiton dispersion, k' = k + p, Akk' = 
Ay (1 + 7p) + A+A_/(l - 7^)1/2, Byy, = A2^7p/(1 + 7p), 
and A-t = 7^"' ± 7k?'' • By averaging first the matrix 
elements A^B in Eq.(5) over the Brillouin zone (they 
are rather regular functions), we evaluate 5Sf^^ ~ 0.15. 
(Another estimation by using the Ising limit for magnons 
gives 0.13). Adding also a conventional 3D spin-wave 
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correction SSgy^ ~ 0.075, one then obtains a staggered 
moment 0.55/xb in fair agreement with experiment. 

The orbital hquid picture offers a simple explanation 
of puzzling Fano-type phonon anomalies observed by Ra- 
man scattering data in insulating titanates (Reedyk et al. 
|p9[). These anomalies are most pronounced in LaTiO^ 
, and we identify their origin as due to the coupling of 
phonons to orbital excitations, noticing that the phonon 
position 300 cm^^) is right in the orbital Stoner con- 
tinuum. 

According to Goodenough-Kanamori rules, a local fer- 
rotype orbital correlations are expected in a spin Neel 
state. A residual interactions (via magnons) between or- 
bitons do indeed produce a triplet pairing of orbitons 
within each fermionic branch. Treating the problem 
on BCS level, we find finite mean-field order parame- 
ter (aj^aj) of p-wave symmetry (specifically, of Px + ipy 
-symmetry for c orbitons, for instance), which opens an 
orbital gap. It is difficult to calculate reliably the value of 
this gap (of the order of a few meV); we predict however 
that a linear 7 -term in specific heat should be released 
when a correlation gap in orbital spectrum is thermally 
washed out. We estimate 7 ^ AOmJ/moleK^. Specific 
heat measurements in LaTiO^ up to T '--^ lOOK would 
be a crucial test for the theory. 

From physical point of view, the proximity of LaTiO^ 
to Mott transition (charge gap is only 0.2 eV is 
actually a key factor, which drives this compound to su- 
perexchange dominated orbital liquid state. This state 
is stable, however, only if the effective orbiton band- 
width, Worb ~ it'^/U, is large enough to suppress JT 
order. When one goes to another end compound, YTiO^ 
(charge gap ^ leV) the bandwidth narrows |16|, thus re- 
ducing superexchange energy scale, and we expect static 
JT ordering of spatially more localized orbitals once the 
ratio Worb /Ejt becomes less than critical. With antifer- 
rotype orbital ordering as observed in YTiO^ a classi- 
cal expectation value of Eq.(2) vanishes, and Hund's cou- 
pling Jh term induces a ferromagnetic state. The compe- 
tition between the frustrated superexchange on one hand, 
and JT plus Hund's couplings on the other hand is ex- 
pected to result in the phase diagram shown in Fig. 3. 

To conclude, we have argued that due to large degen- 
eracy and special geometry of orbitals t2g-superexchange 
system most likely has an orbitally disordered ground 
state. The observed anomalies of LaTiO^ find their nat- 
ural explanations in a proposed orbital liquid picture. In 
metallic manganites, Cg orbital disorder is enforced by 
the presence of mobile carriers . Apparently, the case 
of insulating LaTiO^ is a nice example of that a frus- 
trated t2g orbitals can do the same job alone, without 
doping and well before the Mott transition to a metallic 
state is reached. 
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Prelovsek, and K. Tsutsui are also acknowledged. This 
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FIG. 3. Proposed phase diagram for d}{t2g) Mott insu- 
lator on cubic lattice. At small Hund's and Jahn- Teller cou- 
plings the quantum Neel state with dynamically quenched 
orbital moments is stabilized. Quantum phase transition line 
separates this state from the ferromagnetic phase with static 
orbital/JT order. 
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